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ABSTRACT 


Analytical  relations  have  been  derived  for  calculating  developing  thick, 
axisymmetric,  turbulent  boundary  layers  in  a pressure  gradient  from  two 
simultaneous  differential  equations:  momentum  and  shape  parameter.  An 
entrainment  method  is  used  to  obtain  the  shape  parameter  equation.  Both 
equations  incorporate  the  velocity  similarity  laws  that  provide  a two-parameter 
velocity  profile  general  enough  to  include  any  range  of  Reynolds  numbers. 

Newly  defined  “quadratic”  shape  parameters  which  arise  from  the  geometry 
of  the  thick  axisymmetric  boundary  layer  are  analytically  related  to  the  two- 
dimensional  shape  parameter  by  means  of  these  velocity  similarity  laws. 

The  variation  of  momentum  loss,  boundary-layer  thickness,  local  skin 
friction,  and  local  velocity  profile  may  be  calculated  for  the  axisymmetric 
turbulent  boundary  layers  on  underwater  bodies,  including  the  thick  boundary 
layers  on  the  tails.  The  various  formulations  are  shown  to  correlate  well  with 
available  experimental  data. 

ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  authorized  and  funded  by  the  Independent 
Research  Program  of  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  under 
Project  ZR-023-0101,  Work  Unit  1-1541-002. 

INTRODUCTION 

Accurate  prediction  of  the  viscous  drag  of  bodies  of  revolution  and  the  associated 
boundary-layer  velocity  profiles  requires  solving  equations  for  the  thick,  axisymmetric, 
boundary  layer  on  the  tail.  At  the  high  Reynolds  numbers  of  interest,  the  analysis  involves 
the  usual  difficulties  encountered  with  a turbulent  boundary  layer  in  a pressure  gradient  now 
accentuated  by  the  boundary  layer  being  thick  relative  to  the  body  radius. 

There  are  existing  methods,  having  varying  degrees  of  merit,  for  analyzing  the  thick, 
axisymmetric,  turbulent  boundary  layer  in  a pressure  gradient.  Nelson1  derived  an  energy 
integral  equation  and  used  the  similarity  laws  for  the  velocity  profile.  The  Coles  wake 
function  is  fitted  by  an  awkward  analytical  relation  which  leads  to  complicated  expressions 
for  the  required  boundary-layer  parameters.  The  axisymmetric  energy  dissipation  factor 
which  is  an  essential  element  in  the  solution  is  not  separately  evaluated  but,  instead,  is 


kelson.  D.M.,  “A  Turbulent  Boundary-Layer  Calculation  Method,  Based  on  the  Law  of  the  Wall  and  the 
Law  of  the  Wake,”  U.S.  Naval  Ordnance  Test  Station  (China  Lake,  Calif.),  NavWeps  Report  8510 
(NOTS  TP  3493)  (Nov  1964).  A complete  listing  of  references  is  given  on  page  32. 
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mistakenly  linked  back  to  the  equation  of  motion.  The  result  is  that  the  energy  equation  is 
not  independent  of  the  momentum  equation  which  is  a fundamental  error.  To  correct  this, 
Nelson2  later  empirically  evaluated  the  dissipation  integral  from  two-dimensional  test  data 
whose  application  to  axisymmetric  conditions  has  not  been  established. 

Cebeci  and  Smith3 4  use  a direct  numerical  approach  in  solving  the  partial  differential 
equations  of  motion  for  the  boundary  layer  by  the  method  of  finite  differences.  The 
unknown  turbulent  shear  stress  is  modeled  for  the  most  part  by  an  eddy  viscosity  relation 
previously  obtained  by  Clauser  for  two-dimensional  equilibrium  pressure  gradients.  The  eddy 
viscosity,  originally  related  to  the  two-dimensional  displacement  thickness,  is  extended  by 
Cebeci  and  Smith  to  what  is  essentially  a thin  axisymmetric  displacement  thickness.  At  the 
after  end  of  a body  of  revolution  where  the  radius  of  the  body  becomes  zero,  the  eddy 
viscosity  so  defined  becomes  infinite  and  hence  invalid. 

Patel4,5  developed  a one-parameter  entrainment  method  for  thick  axisymmetric  boundary 
layers  by  extending  the  Head  one-parameter  entrainment  method  for  two-dimensional 
boundary  layers.  In  the  entrainment  equation  derived  by  physical  argument  for  thick, 
axisymmetric  boundary  layers,  Patel  has  assumed  that  the  radius  of  the  body  is  the  length 
scale  required  to  balance  the  equation  dimensionally.  On  the  contrary,  an  analytical  derivation 
later  in  this  paper  shows  that  the  proper  length  factor  is  the  radial  distance  to  the  outside  of 
the  boundary  layer,  which  is  equal  to  the  radius  of  the  body  plus  the  boundary  layer  thickness 
slightly  modified  by  the  slope  of  the  body.  Since  the  radius  of  the  body  becomes  zero  at  the 
very  end  of  the  body,  the  Patel  analysis  shows  no  entrainment,  unlike  the  more  fundamental 
analysis  in  this  paper.  Part  of  the  analysis  of  the  thick,  axisymmetric  boundary  layer  involves 
additional  integrations  of  the  velocity  profile  with  respect  to  the  square  of  the  transverse 
distance  which  may  be  termed  “quadratic”  integration.  For  these.  Patel  uses  approximate 
simple  power-law  velocity  profiles. 


2NeIson,  D.M.,  “Turbulent  Boundary-Layer  Calculations,  Using  a Law  of  the  Wall-Law  of  the  Wake  Method,” 
U.S.  Naval  Ordnance  Test  Station  (China  Lake,  Calif.),  NOTS  TP  4083  (Jun  1966). 

3Cebeci,  T.  and  A.M.O.  Smith,  “A  Finite-Difference  Method  for  Calculating  Compressible  Laminar  and 
Turbulent  Boundary  Layers,”  Journal  of  Basic  Engineering,  Transactions  of  the  American  Society  of 
Mechanical  Engineers,  Vol.  92,  Series  D,  No.  3,  pp.  523-535  (Sep  1970). 

4Patel,  V.C.,  “On  the  Equations  of  a Thick  Axisymmetric  Turbulent  Boundary  Layer,"  Iowa  Institute  of 
Hydraulic  Research  Report  143  (Jan  1973). 

s Patel,  V.C.,  “A  Simple  Integral  Method  for  the  Calculation  of  Thick  Axisymmetric  Turbulent  Boundary 
Layers,”  Aeronautical  Quarterly,  Vol.  25,  part  1,  pp.  47-58  (Feb  1974). 


Gunther6  presents  a complicated  energy  method  that  will  not  be  discussed  herein. 

White  et  at.7  extend  to  thick,  axisymmetric,  turbulent  boundary  layers  a procedure  used 
previously  for  two-dimensional  boundary  layers  which  provides  an  ordinary  differential 
equation  for  the  strcamwise  distribution  of  local  skin  friction  in  lieu  of  the  usual  momentum 
equation.  To  obtain  this  equation,  a constant  shearing  stress  and  the  inner  velocity  similarity 
law  are  assumed  to  hold  across  the  entire  boundary  layer.  These  two  assumptions  are  quite 
contrary  to  experimental  evidence. 

In  this  paper  the  integral  approach  is  followed  wherein  the  axisymmetric  boundary-layer 
equations  are  symbolically  integrated,  and  an  appropriate  velocity  profile  is  inserted  to  obtain 
a numerical  solution.  Integration  of  the  streamwise  equation  of  motion  with  the  assistance  of 
the  equation  of  continuity  gives  the  well-known  momentum  equation.  The  integration  of  the 
equation  of  continuity  gives  a new  entrainment  equation.  The  result  is  an  entrainment  factor 
multiplied  by  the  radial  distance  to  the  outer  edge  of  the  boundary  layer  instead  of  the  radius 
of  the  body  as  used  by  Patel.5 

The  required  velocity  profile  is  supplied  by  the  well-established  velocity  similarity  laws 
for  turbulent  boundary  layers  in  pressure  gradients  which  comprise  an  essentially  two- 
parameter  velocity  profile  system.  A newly  modified  expression  for  the  law  of  the  wake  is 
employed  for  the  analytical  description  of  the  similarity  laws.  The  local  skin  friction  is 
supplied  by  the  similarity  laws. 

The  axisymmetric  boundary  layer  developing  downstream  on  the  tail  of  a body  of 
revolution  becomes  thicker  from  viscous  losses  and  adverse  pressure  gradients  — significantly 
thicker  than  the  body  radius,  which  is  approaching  zero  at  the  tail  end  of  the  body.  What 
happens  is  that  geometrically  the  annular  boundary  layer  in  the  transverse  plane  ends  as  a 
disk  at  the  tail  end.  The  thick  boundary  layer  fills  a thick  annular  space  which  leads  to 
quadratic  terms  in  integrations  of  the  velocity  profile  which  do  not  exist  in  two-dimensional 
flow  or  are  negligible  for  thin,  axisymmetric  boundary  layers.  Newly  defined  quadratic  shape 
parameters  are  accordingly  introduced  and  are  evaluated  from  the  similarity  laws. 

For  simplicity  of  solution,  all  parameters  arc  expressed  in  terms  of  the  usual  two- 
dimensional  momentum  thickness  and  shape  parameter. 

It  is  shown  that  the  momentum  and  entrainment  equations  which  are  ordinary  differen- 
tial equations  may  be  readily  solved  in  a stepwise  fashion  suitable  for  easily  formulated 


6Giinther.  W.,  “Ein  mehrparametrisches  Rechenfahren  fiir  turbulenlc.  inkompressible  Grenzschichten  an 
umstromten  Rotationskorpern,”  Fortschritt-Bcrichte  der  VDI  Zeitschriften,  Series  7,  No.  34  (Oct  1973). 

7White,  F.M.  et  al.,  “Analysis  of  Turbulent  Skin  Friction  in  Thick  Axisymmetric  Boundary  Layers,” 
American  Institute  of  Aeronautics  and  Astronautics  Journal,  Vol.  1 1.  No.  6,  pp.  821  -835  (Jun  1973). 


computer  programs.  Finally,  the  new  formulas  (or  the  quadratic  shape  parameters  are  shown 
to  agree  with  available  experimental  data. 

AXISYMMETRIC  TURBULENT  BOUNDARY  LAYERS 
COORDINATE  SYSTEM 

The  coordinate  system  for  the  boundary  layer  on  a body  o(  revolution  is  obtained  from 
two  families  of  surfaces.  One  family  is  given  by  surfaces  parallel  to  the  given  body  of  revolu- 
tion displaced  by  a normal  distance  v.  The  other  family  is  orthogonal  to  the  family  of 
parallel  bodies  of  revolution  and  is  separated  by  contour  distance  s.  The  coordinate  system 
then  is  s along  a body  meridian  with  s = 0 at  the  nose,  and  y is  normal  to  the  body  surface 
with  y = 0 on  the  body.  For  axisymmetric  flow  only  two  coordinates  are  required. 

At  any  point  (s,  y)  there  is  a radial  distance  to  the  axis  of  the  body  of  revolution  r so 
that  as  seen  in  Figure  1 

r = rw  + y cos  a I ' > 

where  rw  is  the  radial  distance  to  the  body  surface,  and  a is  the  angle  between  the  contour 
tangent  and  the  axis. 


Figure  1 - Thick  Axisymmetric  Boundary  Layer 


EQUATIONS  OF  MOTION 


For  turbulent  flow  the  Reynolds  equations,  which  are  the  averaged  Navier-Stokes 
equations,  are  applied  to  the  boundary-layer  How.  For  the  body-of-revolution  coordinate 
system  the  Reynolds  equations  of  motion®  for  incompressible  flow,  neglecting  longitudinal 
curvature,  are  in  the  s-direction 
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and  in  the  y-direction 
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and  the  equation  of  continuity  is 


0(ru)  0(rv) 

0s  0y 


where  u = velocity  in  s-direction 
v = velocity  in  y-direction 

r = shearing  stress,  including  that  due  to  turbulence 
a,  = turbulent  normal  stress  in  s-direction 
a = turbulent  normal  stress  in  y-direction 
p = pressure 
p = density  of  fluid 


(3) 


(4* 


MOMENTUM  EQUATION 


The  momentum  equation  is  obtained  by  integrating  the  s-equation  of  motion  Equation 
(2)  from  y = 0 to  y = 6 and  incorporating  the  equation  of  continuity  Equation  (4). 

The  result  is  the  momentum  equation  for  thick  axisymmetric  boundary  layers 


d£2  (A*  + 2ni  dU 

ds  + U ds 


- p)  r dy  + 


(5) 


^Granville,  P.S.,  “The  Calculation  of  the  Viscous  Drag  of  Bodies  of  Revolution,"  David  Taylor  Model  Basin 
Report  849  (Jul  1953). 
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12  = momentum  area 


A*  = displacement  area 


n3/S('-o)o,dy 

**■/*(' ■®),dy 


rw  = wall  shearing  stress 

pg  = pressuR'  at  outer  edge  of  boundary  layer 
U = velocity  at  outer  edge  of  boundary  layer 

It  is  to  be  noted  that  no  boundary-layer  simplifications  have  been  made  other  than  neglect 
of  longitudinal  curvature.  The  only  assumption  is  that  a boundary  layer  of  thickness  6 exists 
wherein  all  viscous  effects  are  concentrated. 

ENTRAINMENT  EQUATION 

Besides  the  momentum  equation  Equation  (5),  an  additional  “integral’’  equation  is  needed 
to  calculate  development  of  the  boundary  layer  along  the  body.  The  entrainment  method 
introduced  by  Head9  has  proved  effective  in  this  regard  for  two-dimensional  flows  in  pressure 
gradients. 

The  entrainment  method  will  now  be  developed  for  the  thick  axisymmetric  boundary 
layer.  The  procedure  to  be  followed  is  to  integrate  the  equation  of  continuity  across  the 
boundary  layer  and  then  to  evaluate  the  ensuing  entrainment  velocity  in  terms  of  conditions 
at  the  outer  edge  of  the  boundary  layer.  The  equation  of  continuity  Equation  (4)  is  integrated 
from  y = 0 to  y = 6 


to  give  by  means  of  the  Leibnitz  formula  and  v = 0 at  y = 0 


9 Head,  M.R.,  “Entrainment  in  the  Turbulent  Boundary  Layer,"  Aeronautical  Research  Council  (British), 
RAM  3152  (Sep  1958). 
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d r6  . „ d5 

-J  Ur  dy  - rg  U - + rgVg  = 0 


<l>) 


where  the  subscript  6 refers  to  conditions  at  y = 6.  Now  vg  = v at  y = 6 is  defined  to  be  the 
entrainment  velocity,  and  from  liquation  (1) 


rg  = rw  + 6 cos  a 


(10) 


To  determine  the  entrainment  velocity  vg,  the  equation  of  motion  Equation  (2)  is 
evaluated  at  the  outer  edge  of  the  boundary  layer,  y = 6 


/ /9u\  _ Tbfdp\  l/3rr\  j /dr o,\ 

r§  \3s/g  r&  \ 3y  pUs/5  + P\dy/5  p\bs)t 


(II) 


Various  terms  of  this  expression  are  to  be  examined.  In  general  along  a streamline 


du  _ 3u  + 3u  dy 
ds  ds  9y  ds 


(12) 


and  then  at  y = 6 


/ 9u\  _ dU  / 9u  \ d6 
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Likewise 
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and 


/3p\  _ ^ ^ dU  /9p\  — 

\9s/g  ds  \9y/g  ds 


(15) 


where  the  Bernoulli  theorem  has  been  applied  at  the  outer  edge  of  the  boundary  layer 


(f)  ♦»»?-<> 

\ds/6  ds 

Furthermore  with  r = 0,  and  o,  = 0 at  y = 5 

fdrr)  = (*l)  l^!l\  =n 

\ 9y  /6  Tf>vy)b  \ ds  U 


(16) 


(17) 
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With  the  foregoing  the  entrainment  velocity  given  in  Equation  (11)  becomes 


_ U d6  + _L / dr\  f I faP  + Vl  d6 
U ds  pVWx  pL  3u  Js  ds 


Then  the  entrainment  equation  Equation  (9)  becomes  in  turn 


d r6  r6  (/9r\  p<p  + otr]  d5 

dj0  urdy“7|W6+L'T=_J6* 


The  nondimensional  entrainment  factor  E is  now  defined  as 


_l_(/3r\  , P(p  + qt>]  d6 
pU|\3u/5+L  3u  J5  ds 


This  is  an  expanded  definition  of  E to  include  the  pressure  p variation  and  normal  stress  o, 
variation  in  the  y direction.  Finally  the  entrainment  equation  for  the  thick  axisymmetric 


boundary  layer  is 


if  urdy  = ,{UE 


The  entrainment  equation  Equation  (22)  is  more  fundamentally  based  than  that 
introduced  by  Patel5  for  the  thick  axisymmetric  boundary  layer,  namely. 


urdy  = rwUE 


where  rw  is  used  instead  of  r§.  Patel  adopted  rw  as  a convenient  nondimensionalizing  factor 
in  an  analysis  based  mostly  on  physical  argument. 

For  the  thin  axisymmetric  boundary  layer  5 « rw,  the  entrainment  equation  reduces  to 


dy)’r«UE 


8 


which  is  the  same  as  the  equation  derived  by  Shanebrook  and  Sumner.10 
For  constant  rw  the  entrainment  equation  reduces  finally  to 


d C 8 

-y  udy-Ut 

r* 


Tliis  is  the  same  as  the  Head9  two-dimensional  equation. 

For  the  thick  axisymmetric  boundary  layer  the  entrainment  equation  may  also  be 


expressed  as 


where 


d\J/  dll 

ds  " r§  ' U ds 


’A 


is  a newly  defined  axisymmetric  entrainment  area. 

INTEGRAL  VELOCITY  FACTORS  FOR  AXISYMMETRIC  FLOW 

In  analysis  of  two-dimensional  boundary  layers,  the  displacement  thickness  6*.  the 
momentum  thickness  0 , and  their  ratio  H,  the  shape  parameter,  arc  defined  as 


H = 6 */0  (30) 

For  axisymmetric  boundary  layers,  the  displacement  thickness  5*  becomes  the  displace- 
ment area  A*,  which  has  been  previously  defined  as 


,0Shanebrook,  J.R.  and  W.  J.  Sumner,  “Entrainment  Theory  for  Axisymmetric.  Turbulent,  Incompressible 
Boundary  Layers,”  Journal  of  Hydronautics,  Vol.  4.  No.  4,  pp.  159- 160  (Oct  1970). 


and  the  momentum  thickness  0 becomes  the  momentum  area  SI  which  has  been  previously 
defined  as 


ns/5MHri,y 


Further  analysis  yields  quadratic  factors  such  as  A and  <t>.  From  Equation  (I) 


A*  = rw  6*  + A cos  a 


where 


y dy 


Likewise 


SI  = rw  0 + <t>  cos  a 


where 


**-[  (''n)irydy 


Additional  quadratic  shape  factors  may  now  be  defined  such  as  and  H^, 

A 


Ha=- 


0“ 


and 


<t> 

9 e2 


a 


(• 


c 


c 


Consequently  the  displacement  area  and  the  momentum  area  may  be  related  to 
momentum  thickness  0 by 

A*  = rw  H0  + H^02  cos  a 
and 

SI  = rw  0 + H^02  cos  a 

It  follows  from  the  definitions  given  in  Equations  (I),  (7),  and  (27)  that 


* = rw(«  -«*>  + 


cos  a 


c 


c 


c 


or 


\Jj  = rw  H 0 + 02  cos  a 

where  H is  the  two-dimensional  entrainment  shape  parameter  defined  by 

~ 5-6*  5 

"'TTh 

and  is  a quadratic  entrainment  shape  parameter  newly  defined 


Consequently 


H*  "’(f)3  - + HI2  - H; 


(40) 


(41) 


(42) 


(43) 


SIMILARITY  LAWS  AND  ASSOCIATED  INTEGRALS 


From  a study  of  experimental  data  Nelson*  has  concluded  that  the  similarity  laws  of  the 
turbulent  velocity  profile  apply  to  the  thick,  axisymmetric  boundary  layer  in  a pressure 
gradient  as  well  as  to  the  more  thoroughly  studied  two-dimensional  boundary  layer. 

The  two  velocity  similarity  laws”  implicitly  provide  a two-parameter  system, 
an  inner  law 


and  an  outer  law 


0 < y < y, 


y2  <y,  <y  <6 


(44) 


(45) 


where  uT  = y/rwlp , the  shear  velocity. 

The  inner  law  as  stated  herein  applies  only  to  smooth  surfaces.  For  rough  surfaces  and  for 
drag-reducing  polymer  solutions,  there  are  additional  dimensionless  ratios.  The  outer  law. 
however,  applies  unchanged  to  rough  surfaces  and  to  drag-reducing  polymer  solutions.  The 
two  similarity  laws  overlap  across  the  boundary  layer,  which  leads  to  logarithmic  relations. 11 


’Granville,  P.S.,  “Similarity-Law  Entrainment  Method  for  Turbulent  Boundary  Layers  in  Pressure 
Gradients.”  DTNSRDC  Report  4657  (Dec  1975). 

I I 


Various  new  integral  relations  based  on  the  outer-law  form  are  required  in  the  analysis 
to  come.  The  usual  two-dimensional  relation 


■/'(¥)■< o 


is  now  extended  to  axisymmetric  flow  so  that 


Also  the  usual  two-dimensional  relation 


becomes  now 


■.-rev)'"© 


In  two-dimensional  flow,  the  Rotta-Clauser  shape  parameter  is  defined  as 


where  a = U/uT. 

The  outer  law  may  be  expressed  logarithmically  as 

= - A In  | + B2  (l  - t)*  A q (5,) 

where  is  the  Coles  wake  function,  and  q[|],  the  wake-modification  function,  is  added 

to  make  — = 0 at  y = 6.  A is  a constant  while  B2.  the  velocity-defect  factor,  is  a function 
3y 

* w 

of  s in  arbitrary  pressure  gradients'.  — is  given  by  Moses  as 


! • KI)2  - HI)3 


’^Granville,  P.S.,  “A  Modified  Velocity  Similarity  Law  for  Turbulent  Shear  Flows."  NSRDC  Report  4639 
(May  1975). 


. 
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while  q is  given  by  Granville  as 


(53: 


If  the  outer  law  is  assumed  to  hold  up  to  wall,  0 < y < 6,  which  is  an  excellent 
approximation  at  high  Reynolds  numbers,  the  integral  parameters  become 


B, 


'|  =T2A+T 


4819  213  13  o2 

'2  “ :s:0  A'  * 140  AB2  * 35  B2 


J.  4At3)B2 


and 


443  2 359  4 n 3 „ 2 

J2'^00A‘  + i400  AB2+35B2 


(54] 


(55] 


(56] 


(57] 


where  A is  the  reciprocal  of  the  von  Karman  constant  and  B->  is  the  velocity-defect  factor. 
For  nonequilibrium  pressure  gradients,  B2  = f[s). 


LOCAL  SKIN  FRICTION 


For  the  solution  of  boundary-layer  equations  such  as  the  momentum  equation  Equation 
(5),  a most  essential  ingredient  is  relating  the  coefficient  of  wall  shear  stress  or  local  skin 
I U0 

friction = — to  the  local  Reynolds  number  Ra  = — and  local  shape  parameter  H.  As 

1 12  2 ° v 

pU L o'- 

shown  in  detail  in  Reference  II,  for  two-dimensional  boundary  layers,  this  is  derived  from 
overlapping  the  inner  and  outer  similarity  laws  to  give  in  implicit  form 


0.3462(3.889 

H 


H) 


o + 0.9392  A In  o - A 


In  Rfl  + B| 


2.0938  A 


-A  In  [(H-  l)°-9392/H  1.9392, 


(58 


and  B | is  a law-of-the-wall  (inner)  law  factor  which  is  constant  for  smooth  surfaces  and  a 
variable  for  rough  surfaces  and  drag-reducing  polymer  solutions. 

The  variation  of  a with  H and  required  in  the  next  section  is  derived  in  Reference  1 1 
for  smooth  surfaces  (Bj  = constant)  as 


and 


1 do  2.8885 


a 9H  H(H  - I) 


1 9o 


pl.3462 
L (3. 


(H  - l)o  + AH(H  - 1 
889  - H)a  + 2.7129  A 


.9392)1 

H J 


2.8885  AH 

o3inR0  (3.889  - H)o  + 2.7129  AH 


(60) 


(61) 


ENTRAINMENT  SHAPE  PARAMETER 
AND  ENTRAINMENT  FACTOR 


The  similarity  laws  provide  a two-parameter  evaluation  of  entrainment  shape  parameter 
H,  which  is  derived  in  detail  in  Reference  1 1 as 


Also 


and 


dH  H(H  - 2)  - 


9H  H(H-l) 


H H2  /0. 4739  A 13.7SAz\r  1 , 1 do  1 

~ (H  - 1 ) V,  G + G2-75  / LH(H  + 1 ) ad  Hj 


9H  = / H2  \ /i 

In  Rg  \H  “ 1/  \ 


0.4739  A 13.75  A2 


5 A2\  / 1 do  \ 
.75  J [ a 9 In  R eJ 


(62) 


(64) 


The  entrainment  factor  E,  derived  in  detail  from  equilibrium  pressure  gradients  in 
Reference  1 1 , is 


3 In  R^  \9H/  _ 


AH(H-I)2  1j 

1 1 + P + * j 

_A(H  - l)2  +GH  J ) 

\V  + H P 

J o2  pU2 

where 


G = G + AG 

AG  = - 0.16o  + 3.q2  , 19.7  24.5 

AG  = 0 , o>  24.5 


(66) 

(67) 

(68) 
(69) 
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AG  represents  an  effect  due  to  low  Reynolds  number  which  is  important  in  calculating  the 
turbulent  boundary  layer  just  after  transition  from  laminar  flow. 


QUADRATIC  SHAPE  PARAMETERS 


From  previous  definitions,  the  quadratic  displacement  shape  parameter  may  be 
written  in  the  form 


(0/6)2 


where 


HV  §)(*)-&) 

and 

H'KK-SMi) 

Evaluation  of  depends  on  the  choice  of  velocity  profile  First  a simple  power-law 

velocity  profile  will  be  used  as  a basis  of  comparison.  Then  a two-parameter  velocity  profile 
arising  from  the  velocity  similarity  laws  will  be  utilized  to  derive  a relation  for  H^. 

The  power-law  velocity  profile  is  given  by 

Ht)" 

where  n is  a constant. 

Elementary  integration  in  Equations  (71)  and  (72)  results  in 

_ h2(H+  n2 

2(H  - I)  (H  + 3) 

where  H obtained  from  integrations  of  Equations  (28)  and  (29)  is  substituted  for  n. 

From  the  velocity-similarity  laws.  is  derived  as  follows  from  the  definition  of  J,  in 
Equation  (47)  and  from  Equation  (71) 

a J. 


I >• 


15 


and  - is  given  by  Equation  (72). 
o 

Use  of  the  power-law  velocity  Equation  (73)  results  in 


H - h2<h  + , o , 

* (H  + 3)  (H  - I)  H + I 

For  the  velocity  similarity  laws  is  derived  as  follows.  From  the  definition  of  Jj  in 
Equation  (47)  and  J,  in  Equation  (49) 

0 J1  J2 

t (X! 

5“  0 a2 

0 

Use  of  t in  Equation  (76)  and  Equation  (78)  for  o results  in 
o 


«OT&,&)-h20) 


— — has  already  been  evaluated  in  E()uation  (80). 


— r is  related  to  G by  means  of  Equation  (57)  and  (54)  with  B->  as  the  implicit  parameter. 

'■  , 

J2  12  J 2 

At  the  limiting  condition  of  separation  B-,  -*■  and  —z.  rjr.  -7  as  a function  of  G gives 

1 2 35  | - 

. > I 

the  approximation 


J2  12  0.1821  A 


2 35 


which  satisfies  the  condition  at  separation  where  G -*■  Using  the  relationship  between  0 
and  G in  Equation  (78)  results  in 


0.1028  H2  (H  + 3.336)  + /p.  1 821  AH3\  1 + f 3 A H6  1 J_ 
H - I \ H-l  o Lh  - n4l  „3 


The  quadratic  entrainment  shape  parameter  is  given  by  Equation  (43)  as 


H - ' 6 H 
H0  " T 0 
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For  power-law  velocity  profiles  Equation  (73) 


6 _ H(H  + 1) 
9 H - 1 


Then  from  Equation  (74)  for  H^,  it  follows  that 


2H2(H+I)2 


(H  - I)2  (H  + 3) 


For  the  velocity  similarity  laws 


0474A  + iAi 

G r.2.7: 


using  Equations  (41)  and  (62). 
Finally 


3 AH6  _1_ 
(H-  l)4  o3 


from  Equations  (78),  (81),  (43),  and  (91).  ' 

The  quadratic  shape  parameters  H^,  H^,  and  are  shown  in  Figures  2 through  4. 
Calculations  of  the  similarity  law  formulations  are  based  on  values  of  the  constants  for  the 
similarity  laws  as  given  in  Reference  I 1.  Here  the  values  of  the  constants  are  given  by 
A = 2.606  and  B,  = 3.88,  which  are  chosen  so  that  the  drag  coefficient  closely  fits  the  well- 
accepted  Schoenherr  formula  for  smooth  fiat  plates  for  the  case  of  zero  pressure  gradient. 
Since  the  similarity  laws  provide  a two-parameter  variation,  the  quadratic  shape  parameters 
vary  with  Reynolds  numbers  as  well  as  H.  The  plots  show  the  variation  of  H^,  H^,  and 
at  two  momentum  thickness  Reynolds  numbers,  R^  = 103  and  105.  It  is  to  be  noted  that 
the  greatest  effect  of  Reynolds  number  variations  is  at  the  lower  values  of  H.  For  comparison 
the  power-law  relations  are  also  plotted.  The  agreement  at  high  Reynolds  numbers  is 
surprisingly  close. 


Experimental  values  obtained  by  Patel  et  al.13’14  are  also  plotted.  Values  of  Ha,  Ha, 


and  H^  for  a thin,  axisymmetric  boundary  layer  — < 0.2  are  not  plotted  inasmuch  as  the 


Patel,  V.C.  et  al.,  “An  Experimental  Study  of  the  Thick  Turbulent  Boundary  Layer  Near  the  Tail  of  a Body 
of  Revolution,”  Iowa  Institute  of  Hydraulic  Research  Report  142  (Jan  1973). 


Patel,  V.C.,  “Measurements  in  the  Thick  Axisymmetric  Turbulent  Boundary  Layer  Near  the  Tail  of  a Body 
of  Revolution/*  Journal  of  Fluid  Mechanics,  Vol.  63,  part  3,  pp.  345  -367  (1974). 
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SHAPE  PARAMETER  H 


SHAPE  PARAMETER  H 


quadratic  effect  is  minor.  The  values  of  HA,  H^,  and  at  the  high  val  'e  of  H = 2.23  are 

suspect  because  at  the  99-percent  station,  the  boundary  layer  is  undergoing  transition  to  a 

wake.  In  general,  it  should  be  noted  that  the  agreement  between  the  experimental  values 

and  the  predicted  values  is  excellent  for  and  H^. 

In  the  case  of  the  experimental  values  are  uncertain  because  of  the  functional 

dependence  of  H^,  on  the  values  of  boundary-layer  thickness  5.  In  general,  the  evaluation  of 

boundary-layer  thickness  is  an  imprecise  process  unlike  the  evaluation  of  boundary-layer 

momentum  thickness  0 which  involves  an  integration.  To  illustrate  this,  the  values  of  - are 

6 0 

examined  in  Figure  5.  Also  given  are  values  of  - obtained  from  the  Head  formulation  of  the 

~n  0 

entrainment  shape  parameter  H 1 1 


H = 1.535  (H  - 0.7r2-715  + 3.3  (93) 

and 

|=H  + H (94) 

from  Equation  (41). 

As  seen  in  Figure  5,  the  experimental  values  of  Patel  et  al. 13,14  are  uniformly  higher  than 
the  computed  values  of  -.  This  may  be  due  to  the  uncertainly  in  determining  the  value  of 
boundary-layer  thickness  which  is  not  a sharply  defined  quantity. 


METHOD  OF  SOLUTION 

For  a given  geometry  rw|s]  and  external  flow  velocity  U|s|,  development  of  the 
boundary  layer  may  be  calculated  from  a stepwise  solution  of  the  momentum  equation  and 
the  entrainment  equation.  The  primary  dependent  variables  to  be  calculated  are  0 and  H. 

All  other  factors  are  functions  of  9 and  H.  The  momentum  area,  displacement  area,  and 
entrainment  area  are 

ft  = rw  0 + H^fl2  cos  a (37) 

A*  = rw  H0  + H^02  cos  a (38) 

i//  = rw  H0  + H^02  cos  a (40) 

Differentiating  f2  and  \p  with  respect  to  9 and  H gives 
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Figure  5 - Ratio  of  Boundary -Layer  to  Momentum  Thicknesses 


(93) 


['« ♦,(co,«"(2 H* * Fi^)]  M * ['C05“,#2  7ff]dH 


, ->  dMd>  dU 

= dft-0drw  - HqO2  d(cosa)-(cosar)0‘  3^^  JJ 


3 Ha  3 Ha 

dH*=  Th  dH  * mrfe  dinR» 


dO  dU 

d ln  = J + “JJ 


[r»  (H  * Fhf^)  ♦•(«»«»  "*(2  ♦ 5^F9)] de  * [r* 6 IS 4 ,cos 92  TFr] d h 

2 / d H \ 2 

= d^-H0drw-rwF^  ^dU-H^02d(cosa)-(coso()^g-i^j^dU  (96) 


dfi  and  d <//  are  obtained  from  the  momentum  equation  and  entrainment  equation 

dn  = (r"^)ds‘!A^dU  (,7) 

and 

d*^  = [rw  +0)(?  cosa]Eds--^-dU  (98) 

Then  Equations  (93)  and  (96)  are  solved  as  simultaneous  algebraic  equations  for  dd  and  dH 
for  each  differential  step  ds. 

The  momentum  equation  Equation  (5)  is  solved  without  including  the  transverse 
pressure  and  normal  stress  terms  inasmuch  as  reliable  formulations  do  not  exist.  Patel5 
recommends  the  use  of  the  measured  pressure  distribution  on  the  surface  of  the  body  instead 
of  pg  in  determining  U in  Equation  (16)  as  compensation  for  excluding  the  transverse 
pressure  term.  The  normal  stress  term  is  usually  negligible,  except  close  to  separation. 


SUMMARY 


A new  entrainment  equation  is  derived  for  thick  axisymmetric  boundary  layers 


New  quadratic  shape  parameters.  HA.  H0.  and  are  introduced  and  evaluated  by  the 
velocity-similarity  laws  in  two-parameter  relationships,  namely,  as  functions  of  H and  R0. 
Agreement  with  test  data  is  excellent.  This  tends  to  corroborate  the  assumption  that  velocity 
similarity  laws  apply  also  to  thick,  axisymmetric.  turbulent  boundary  layers. 

All  the  necessary  analytical  relations  have  been  derived  for  a two-parameter  calculation 
of  thick,  turbulent  boundary  layers  in  pressure  gradients.  The  momentum  loss,  the  boundary- 
layer  thickness,  and  the  velocity  profile  may  be  determined  at  each  station  in  the  downstream 
direction.  All  the  pertinent  relations  are  assembled  in  the  appendix  for  numerical  calculation 
of  the  thick,  axisymmetric.  turbulent  boundary  layer  and  the  associated  velocity  profiles. 
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APPENDIX  A 

CALCULATION  PROCEDURE 


All  the  necessary  analytical  relations  for  calculating  development  of  the  thick,  axisym 
metric,  turbulent  boundary  layer  and  the  associated  velocity  profiles  are  now  listed. 

The  pertinent  equations  have  been  nondimensionalized: 


L = length  of  body 
G = axial  distance  from  nose 


GIVEN 


SOLUTION 


solve  as  simultaneous  equations  for  6 


and  6H  where  6 is  a differential  increment 


•r*mm  i*f>  rm 


3H0  0.20S6H(H2  + 0.168  H-  3.336)  AH3  [b.1821  ( 9H3  ~[  / 1 

3H  " (H-l)2  ~ (H  - I ) L o (H-  I)3  o3J  3H/ 

0.1821  AH2(2H  - 3)  + 6AHS(H  - 3) 

(H  - 1 )2  o (H-1)5o3 

3H0  _ AH3  [“0. 1821  , 9 H3  1 A do  \ 

3 1nRe  "'(H-  l)|_  o %H-  l)3  o3J  V°  3ln,V 

A = 2.606 
B,  = 3.88 

aj  - 6.1,  jj  ~ 1-^1.  ^3  = 1.37 


INITIAL  CONDITIONS  AT  TRANSITION  (Reference  11) 


/ uTy\  uTy  ur6 

Velocity  profile  for  ( ~ ) < — < — 


L /L  v v 


U / UT^  ^ W 


7=3(S)2-4)3 


=(s)2M) 


J = B | + AlnA-A 


— ) = B,  + A In  A 


Velocity  profile  for  0 < 


s*  <(ujA 

v \ ^ 


U UT^ 

laminar  sublayer  — = — . 

uT  v 
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